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ABSTRACT 

We revisit the well known discrepancy between the observed number of Milky Way (MW) dwarf 
satellite companions and the predicted population of cold dark matter (CDM) sub-halos, in light of 
the dozen new low luminosity satellites found in imaging data from the Sloan Digital Sky Survey 
(SDSS) and our recent calibration of the SDSS satellite detection efficiency, which implies a total 
satellite population far larger than these dozen discoveries. We combine a detailed dynamical model 
for the CDM sub-halo population with simple, physically motivated prescriptions for assigning a 
stellar content to each sub-halo, then apply observational selection effects and compare to the current 
observational census. Reconciling the observed satellite population with CDM predictions still requires 
strong mass-dependent suppression of star formation in low mass sub-halos: models in which the stellar 
mass is a constant fraction F»(f2(,/J7 m ) of the sub- halo mass M sat at the time it becomes a satellite fail 
for any choice of F*. However, previously advocated models that invoke suppression of gas accretion 
after reionization in halos with circular velocity V c j rc < V^ r it ~ 35 km s _ can reproduce the observed 
satellite counts for —15 < My < 0. Successful models require F* « 1CP 3 in halos with V c - uc > V" C rit and 
strong suppression of star formation before reionization in halos with T4; rc < 10 km s" 1 ; models without 
pre-reionization suppression predict far too many satellites with — 5 < My < 0. In this successful 
model, the dominant fraction of stars formed after reionization at all luminosities. Models that match 
the satellite luminosity distribution also match the observed heliocentric radius distribution, and they 
reproduce the observed characteristic stellar velocity dispersion cr* ~ 5 — 10 km s -1 of the SDSS 
dwarfs given the observed sizes (~ 50 — 200 pc) of their stellar distributions. The model satellites 
have M(< 300 pc) ~ 10 7 M Q as observed even though their present day total halo masses span more 
than two orders of magnitude; the constancy of central masses mainly reflects the profiles of CDM 
halos. Our modeling shows that natural physical mechanisms acting within the CDM framework can 
quantitatively explain the properties of the MW satellite population as it is presently known, thus 
providing a convincing solution to the 'missing satellite' problem. 

Subject headings: Galaxy: halo - Galaxy: structure - Galaxy: formation - Local Group 



1. INTRODUCTION 

The inflationary cold dark matter scenario predicts an 
initial fluctuation spectrum with power that continues 
down to small scales, and in consequence it predicts a 
mass function of dark matter halos that rises steeply to- 
wards low masses. A significant fraction of these ha- 
los survive as gravitationally self-bound units long after 
falling into more massive h alos . As pointed out f orce- 
fully bv lKlvpin et~aTl (|1999l ) and lMoore et alJ (|1999l ). the 
predicted number of sub-halos within a Milky Way-like 
galaxy halo greatly exceeded the then known numbers of 
Milky Way or Local Group dwarf satellites, when sub- 
halos and observed dwarfs were matched based on ve- 
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locit y dispersion or corresp onding circular velocity (see 
also iKauffmann et al.l ll993). This discrepancy between 
predicted and observed numbers has become known as 
the "missing satellite problem." 

Proposed solutions fall into three general cate- 
gories. The first modifies the properties of dark 
matter or the primordial fluctuations from inflation 
in a way that eliminates th e low mass dark mat - 
ter s u b-halos themselves (e . g . IKamionkowski fc Liddl 
20051 iSpergel fc Steinhardt] 120001 : iBode et all 120011 
Zentner fc Bullock! 12003 )7" The second appeals to as- 



trophysical mechanisms that suppress star formation in 
low mass halos so that they do not become observ- 
able dwarf satellites; photo-heating by the meta-galactic 
UV background is an attractive mechanism because it 
naturally introduc es a cutoff at approximately the cor- 
rect velocity scale (|Bullock et al J 120001 : iSomer vm~[2002; 
iKravtsov et al.ll2004f ). The third possibility, arguably a 
variant of the second, is that the numerous dwarf com- 
panions of the Milky Way actually exist but have been 
missed by observational searches. 

In this paper we revisit the "missing satellite prob- 
lem" with particular emphasis on the role of the new 
dwarf companions discovered in im aging data from 
the Sloan Digital Sky Survey (SDSS; lYork et al.ll2~)0~ : 
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lAdelman-McCarthv et all 120081) . There are now about 
a doz e n of t he se (iWillman et al.l 120051: IBelokurov et alj 



2001. [20071: [Zucker et all 120061: llrwin et all I2007f: 



Koposov et al.1 



2007H Walsh et all 120071 : a couple of sys- 



terns still have ambiguous status), most of them at least 
an order of magnitude less luminous than the faintest 
of the previously kno wn, "classical" sate l lites. 7 Spectro- 
scopic follow-up (e.g. IMartin et al.ll2007l : iSimon k. Gehal 
2007; IGeha et al.l [2008) for many of them indicates that 
they are indeed dark matter dominated systems, even 
though most are fainter than ty pical globular clusters , 
as low as only ~ 1000 L Q (e.g.. IBelokurov et all 120071 : 
IMartin eT al. 2008). Remarkably, almost all of the newly 
found faint satellite galaxies have stellar velocity dis- 
persions in the range 3 — 10 km s~ , though their lumi- 
nosities vary widely. Similarly, the total masses within 
the inner 300 pc sp an less than an order of magnitude 
(IStrigari et al.N2008l) . 

Since the SDSS imaging in which these satellites have 
been discovered covers only ~ 20% of the sky, a naive ac- 
counting would increase the estimated number of Milky 
Way companions by 5 x 12 = 60, in addition to the ten 
classical satellites. However, IKoposov et al.1 (|2008f) use 
a well-defined identification algorithm to show that the 
SDSS dwarfs are also subject to strong radial selection 
effects. Most of the newly discovered objects could only 
have been found within distances of 50-100 kpc, much 
smaller than the inferred virial radius of the M ilky Way's 
dark matter halo (~ 280 kpc for p V i T /p — 340: lXue et al.l 
120081) . The faintest SDSS dwarfs are detectable over 
only 1/1000 of the halo virial volume (inclu d ing th e fac- 
tor of five for sky coverage ). IWalsh et al.l (|200l) have 
recently reached similar conclusions based on an inde- 
pendent identification algorithm and independent Monte 
Carlo tests. 

Such analyses are the basis for 'volume corrections' for 
the faint Milky Way satellite population. With proper 
volume corrections applied, the luminosity function of 
faint Milky Way satellite galaxies turns out to be a 
rather shallow power law in the range —15 <My < —3 
(|Koposov et al.ll20"08l ). These results in turn imply that 
the number of satellites brighter than My = — 3 is ~ 
80 or more, an d the number above My = could be 
a few hundred. I Toller ud et al.l (|2008D reached a similar 
conclusion, adopting a radial sat ellite distribution base d 
on the Via Lactea simulation of IDiemand et a l. ( 2007). 
Even this census counts only dwarfs that are above the 
effective surface brightness thre shold for SDSS detection. 
With the lKoposov et all (|2008l ) detection algorithm, this 
threshold is approximately 30 mag arcsec -2 (T^-band), 
with a weak dependence on luminosity and distance. 
The dwarfs found in SDSS have surface brightnesses that 
range from 24 to 30 mag arcsec -2 . 

Studies of the high redshift Lya forest indicate that 
the small scale power expected in the standard ACDM 
scenario (inflationary cold dark matter with a cosmolog- 
ical constant) is indeed present in the primordial fluctu- 
ation spectrum dNaravanan et al.ll2000l : iViel et all 120051 : 
I Abazajianl 120061 : ISeliak et al.l 120061 ). Astrophysical sup- 
pression of star formation, and photo-ionization suppres- 
sion in particular, has emerged as the most plausible and 



hence popular solution to the "missing satellite" conun- 
drum. Within this category, there have been different 
proposals about what sub-halos host the observed dwarf 
satellites. IBullock et al.l ((200 03 suggested that the ob- 
served dwarfs are those whose sub-halos assembled a sub- 
stantial fraction of their mass before reionization, and 
thus before t he on set of photo-ionization suppression. 
IStoehr et al.l (|2002fl suggested that the measured stel- 
lar velocity dispersions are well below the virial velocity 
dispersions of the dark matter sub-halos, and that the ob- 
served dwarfs occupy sub-halos that are still above the 
veloc ity threshold where st ar-formation suppression oc- 
curs. iKravtsov et alj (|2004l ) used N-body simulations to 
show that roughly 10% of sub-halos lose a large fraction 
(~ 90%) of their mass during dynamical evolution with- 
out being completely disrupted; they suggested that the 
observed dwarfs occupy sub-halos that were above the 
suppression threshold at the time they became satellites 
but have suffered extens ive mass loss since then. These 
pape r s and others (e.g., iSomervilld [20021 : iStrigari et al.l 
120071 : lOrban et alj I2008D focus on explaining the clas- 
sical (pre-SDSS) dwarf spheroidal population, with lu- 
minosities in the range —8 < My < —15 (excluding 
the Magellanic clouds) and stellar velocity dispersions 
in the range 8 km s~ < cr* < 25 km s~ . The recently 
discovered SDSS dwarfs have much lower luminosities 
(—8 < My < —1.5), lower surface brightness, and some- 
what lower velocity dispersion (o% ~ 5 km s~ ), so they 
could have a distinct formation mechanism, or they could 
form a continuum with the classical dwarf sphcroidals. 

The new SDSS discoveries and their quantified de- 
tectability are the basis for the model-data comparison 
in this paper. We construct and test models of the 
Milky Way dwarf satellite population that incorporate 
Monte Carlo realizations of merger trees for 10 12 Af Q 
(main galaxy) halos, a detailed analytic model for the 
dynamical evolution and disruption of sub-halos, and a 
variety of recipes for assigning stellar masses to these sub- 
halos motivated by ideas in the existing literature. For 
most of our models, we assume that a sub-halo can only 
accrete gas to form stars (a) before the epoch of reioniza- 
tion or (b) after reionization if its virial velocity exceeds 
a critical threshold before it enters the Milky Way halo 
and becomes a satellite. The s pirit o f the exercise is sim- 
ilar to that of IBullock et al.l (|2000T ). but the dynamical 
modeling of sub-halos is more sophisticated, and we are 
now in a position to include directly the (strong) con- 
straints imposed by the SDSS dw arfs accounting for th e 
radial selection function found bv lKoposov et all (|2008f) . 
In contrast to most previous studies, we treat the lumi- 
nosity distribution as the primary test of models, rather 
than the stellar vel o city dispersions or central masses 
(IStrigari et alj [20071 12001 ILi et all 12001 IMaccio et all 
2009) , or the inferred but unobservable sub- halo circular 
velocities. This emphasis is motivated by the fact that 
the luminosity is the foremost quantity that matters for 
the observational selection. We consider stellar velocity 
dispersions and central masses as an additional test, but 
their interpretations are affected by the uncertainty in 
the dark matter profiles of the sub-halos associated with 
observed dwarfs. 



7 Throughout the paper we use "faint" and "bright" to refer to 
intrinsic luminosity rather than apparent brightness. 



2. THE POPULATION OF DARK MATTER SUB-HALOS IN 
THE MILKY WAY 



Milky Way satellite galaxies. 
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Our model for the Milky Way satellites is based on 
the cold dark matter scenario, with each satellite form- 
ing initially in a separate dark matter halo that at some 
point falls into the Milky Way's dark matter halo. We 
refer to the bound dark matter satellites orbiting in the 
Milky Way halo as sub-halos. A sub-halo may or may 
not correspond to a dwarf satellite galaxy, depending on 
whether it contains an observable number of stars. In 
this Section we describe our model for computing the 
dynamical evolution of sub-halos. 

We use the dyna mical dark - matte r-only model of sub- 
halos developed bv lYoo et alJ (|2007l ) to compute the sub- 
halo population and its orbital dis t ribut ion. This model 
is described in detail in lYoo et al.l (|2007| ). where a much 
larger halo of 10 15 M was considered as a model of a 
massive cluster of galaxies. Here we consider instead a 
final halo of 10 12 M Q at the present time as a representa- 
tion of the Milky Way galaxy. Despite the change in the 
final halo m ass, the mo d el rem ains basically the same as 
described in lYoo et all (|2007h , so here we make only a 
brief summary of its description. 

The model uses the extended Press-Schechter formal- 
ism to generate a Monte Carl o merger tree of the par- 
ent halo at the present time (|Press fc Schechterl [l97l 
iBond et "al~lll991f ). We follow the dynamical evolution of 
all the sub-halos with masses Mh > 10 6 M Q until they 
merge with the Milky Way and lose their mass below 
Mh = 1O 5 M . All halos start as isolated objects, and 
they grow in mass by accretion and mergers for as long as 
they remain isolated. At some redshift, z sa t, they merge 
into a larger halo (either the Milky Way or another ob- 
ject that will become a Milky Way sub-halo). After this 
merger, the object has become a satellite or sub-halo 
and it stops growing in mass. It can subsequently lose 
mass by tidal stripping when it passes near the center 
of its parent halo or undergoes encounters with other 
sub-halos. The sub-halo is subject to dynamical friction, 
which tends to shrink its orbit, and to random encounters 
with other sub-halos, which on average expand the orbit. 
The orbital eccentricity is also subject to random vari- 
ations. The model allows for the presence of sub-halos 
within other sub-halos. When a sub-halo is disrupted, 
any sub-halos it contained are dispersed into the new, 
larger parent halo. This simple analytic model is able to 
reproduce the sub-halo mass function, in reasonably good 
agreem ent with that found in numerical iV-body simu- 
lation s (jZentner et al.ll2005l : IShaw et al.l [20061 : lYoo et all 
l2007f l. For the present purpose, this approach has the 
advantage (over N-body) of easily affording the required 
mass resolution and multiple halo realizations. 

We adopt a flat ACDM cosmology with matter density 
fi TO = 0.24, baryon density fib = 0.04, power spectrum 
normalization erg = 0.8, Hubble constant h = 0.7, and a 
primordial spectral i ndex n K = 0.95, consistent with re- 
cent measurements (jSpergel et al.1 [20071 iTegmark et al.1 
2006). The matter power spectrum is computed b y us- 
ing the transfer function of Eise nstein fc Hu (1999J). We 
generate six Monte Carlo merger trees of a Milky- Way 
sized halo. Each realization provides the sub-halo mass 
function, their orbital elements and density profiles at 
the present time. Our statistical results are the average 
of the six different realizations. 

The dynamical model of Yoo et al. (2007) uses the Jaffe 



profile and its velocity dispersion to model sub-halos and 
their dynamical interactions, for reasons of numerical 
simplicity and because large galaxies that are tidally- 
limited satellites of a larger halo are reasonably well mod- 
elled by a Jaffe sphere for their baryon plus dark matter 
density profiles. However, the very low-mass dwarf satel- 
lites tend to be dominated by dark matter even in their 
inner parts. We therefore make an adjustment to better 
connect our Monte Carlo simulation results to the ob- 
served Milky Way dwarf galaxies: we use the sub- halo 
masses and orbital elements, which are the quantities 
most robustly computed in the Yoo et al. (2007) model, 
but we calculate the density profiles and velocity dis- 
persion s of sub-halos assum ing that they have an NFW 
profile (|Navarro et al.| [l997V Using the standard spher- 
ical collapse model, the virial radius of an isolated halo 
is assigned as 



7?v 



3M 



halo 



4-7T A c p m (z) 



(1) 



where A r 



(187T 2 



82 x 



39ai 2 )/(l + x) 

(|Brvan fc Normanl Il998h . x = -(1 - fl m )/(Q m (l + 
z) 3 + 1 — fl m ), and the mean cosmic density is 
Pm(z) = fl m p c (l + z) . The halo concentration c is 



Bullock et al. (2001), 



computed using the relation from 

scaled to as = 0.8 according to iMaccio et al.l (|2007i ). 

with c = 0.8 x 9 x (AVIO^/i^Mq)" ' 13 /(I + z). 

For the model in this paper, we use in particular the 
sub-halo masses at two different special epochs: M re i = 
-AT t ot(z = z rc j) when the universe reionizes and the photo- 
ionization background starts to suppress the star forma- 
tion efficiency in low mass halos, and M sat = M tot (z sa t) 
at the epoch when a halo merges into a larger halo and 
we presume that subsequent star formation and gas ac- 
cretion is halted in the sub-halo. We shall also use below 
the halo circular velocity T4; rc , which is the virial circu- 
lar velocity, V chc = [(G M tot )/i? v ir] 1/2 - (Here M tot refers 
to the total mass including dark matter and a universal 
fraction of baryons.) 

In Figure [1] we show the distribution of M su b-haio at 
redshifts z = (left panel), z = z sat (middle panel) and 
z = 8, II, and 14 (right panel). As expected, the mass 
distribution is close to a power-law, except near the res- 
olution limit of our simulations. Also, in Figure [2] we 
show the accretion history of the MW sub-halo popu- 
lation, by plotting the halo masses of the present day 
MW sub-halos at the time of accretion vs. the red- 
shift at which they became satellites of larger halos. We 
see that most of the MW sub-halos became satellites at 
z < 2. Most of the accreted satellites have small circu- 
lar velocities V c i IC < 20 km s~ , so they lie in a range 
where gas accretion and star formation ar e likely to be 
supp r essed after the epoch of reionization (lOuinn et aTl 
119961: iThoul fc Weinberg! 1 19961: iBullock et al.H200Cf) . 

3. POPULATING THE DM HALOS WITH STARS 

3.1. Recipes to assign stellar masses to sub-halos 

To make direct observational predictions from these 
models, we populate each sub-halo in a given Monte- 
Carlo realization with stars according to a sequence of 
recipes, then test how many of these satellites could have 
been found within the SDSS. Some of these recipes are 
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Fig. 1. — Mass distribution of dark matter sub-halos at different epochs. Left panel: The present-day {z = 0) mass function of sub-halos 
within a Milky Way-like halo. This histogram is repeated in the other two panels for reference. Middle panel: Distribution of mass that 
present-day sub-halos had at z = z sa t, the epoch at which they became a satellite within a larger halo {thin line); tidal stripping of satellite 
halos is manifesting important. Right panel: Mass distribution of present day MW sub-halos at the epoch of reionization, for z rc i = 8 
{dotted), 11 {thin solid), and 14 {dashed). All panels reflect the average of six different realizations of MW-like halos. The flattening below 
At = 10 6 /i _1 M Q and the sharp cut-off at At = 10 s h' 1 M© arise from the mass resolution limits of our simulations. 




Fig. 2. — Epochs when (sub-)halos were accreted into larger ha- 
los, and masses at that time. This Figure illustrates the results 
from one Monte-Carlo realization of the semi-analytic model, with 
each point showing the redshift z S at at which a sub-halo first be- 
came a satellite in a larger halo against its total mass M sa t at 
that epoch. The small panel on the right shows the distribu- 
tion of 2 S at ■ Solid and dashed lines show the locus of halos with 
K:irc( 2 sat) = 40 km s — 1 and 20 km s — , respectively. 

mathematically simple illustrations, while others are mo- 
tivated by the expected effects of ionization and cooling 
physics as discussed in the introduction. For reference, 
the nomenclature of the recipes is summarized in Ta- 
ble [1] In all cases we calculate the stellar mass based 
on the sub-halo mass (dark matter plus baryons in the 
universal fraction) at the accretion epoch z sat , denoted 
-Msat- We implicitly assume that satellites do not ac- 
crete new material to form additional stars and that tidal 
stripping of the dark matter does not affect the stellar 



content of the satellite if it survives to the present day. 
Simulations suggest that these assu mptions are reason- 
able but not perfect ap proximations (|Simha et alj 120081 : 
iPenarrubia et al.ll2008fh 

We begin with the simplest model (denoted Model 1 A) , 
that the stellar mass is a constant fraction of the sub- halo 
mass at the time of accretion into the main halo: 

M„ = /* x M sat . (2) 

The a rguments of iKlvpin et al.l (|1999l ) and iMoore et al.1 
(1999) suggest that this model will fail badly, and we 
show that it does indeed fail despite the new satellite dis- 
coveries and the radial selection biases that affect them. 
There is ample evidence that the efficiency of star forma- 
tion declines rapidly towards low m asses even well above 
the d warf satellite regime (e.g., Ivan den Bosch et alJ 
[2007h . In Model IB, we allow the stellar fraction to vary 
as a power law of M sat below a threshold Mq: 

M. = /* X min ( (^j^J , A x M sat . (3) 

Our second approach to modeling stellar masses in- 
cludes the effects of a pervasive energetic radiation 
field after the epoch of reionization, which heats gas 
and hence keeps it from accumulating at the cen- 
ters o f low -mass halos. Calc u lation s by iQuinn et alJ 
(|1996t ) and iThoul fc Weinberg! (|1996f l showed that gas 
accretion in halos with the circular velocities below 
Vciic ~ 30 — 40 km s _1 is strongly suppressed, while sub- 
stantially larger halos are m i nimal l y affected (see also 
iWeinberg. Hernquist. fc Katzl Il997t iGnedinl |2000T) . In 
this spirit, we assume that halos below a critical circu- 
lar velocity form no stars after reionization, and we thus 
assign stellar masses 

M = V* X Msat if Vcirc ^ Zsat ^ > Vcrit (4) 
1 /* X M rei if V ciTC ( ) < Vcrit- ' ' 

This model (Model 2) has three adjustable parameters 
~~ f*i Vcrit, and z ro ; — with expectations that Vcrit ~ 
20 -40 km s" 1 and z rei ~ 11 fe.g. lWeinmann et al.ll2007l ; 
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Dunkl ev et alj|200£ ) . The approach is similar to that of 



Bullock et al.l (|2000f ). except that we treat V cr n as free pa- 



rameter, and the stellar mass formed before the epoch of 
reionization is assigned using M* — /* x M re i, instead of 
simply dividing galaxies into "observable" or "unobserv- 
able" classes based on the fraction of the mass accreted 
by Zrci- 

Our third class of models is similar to the second, but 
it replaces the sharp threshold of equation ((4]) with the 
continu o us tra n sition found i n num erical simulations by 
IGnedinl (f2000h . IHoeft et alj (pOOtl) , and IQkamoto et~afl 
(2008). The numerical results in these papers ca n be de- 
scribe d fairly well by a formula similar to that in IGnedinl 
(2000), with the fraction of baryons that cool in low mass 
halos suppressed by a factor [1 + 0.26(V r cr i t /V^i rc ) 3 ]~ 3 ; 
well after the reionization redshift, the critical velocity 
is found to be approximately independent of redshift. 
IGnedinl (|2000D found y crit ~ 40 km s~\ but these re- 
sults were artificially affected by n umerical resolution 
(N. Gnedin, private communication). IHoeft et"aLl (|2006l ) 
and IQkamoto et alj (|2008l ) find V clit ~ 25 - 30 km s" 1 . 
Including the pre-reionization contribution to M*, this 
model (Model 3A) becomes 



M* = 



/. x (M sat - M Tci ) 

(l+0.26(T/ CI ,t/Kirc(2 S at)) 3 ) : 



+ /, x M rei . (5) 



The assumption that all halos can form stars before 
z re i may not be justified because in halos with virial 
temperature T vir < 10 4 K (Kirc ^ 10 km s -1 ) the gas 
does not get hot enough to cool by atomic processes, 
and simulations that include molecular cooling suggest 
that gas cooling and star formation i s very inefficient 
in su c h halos (IHaiman et alj Il997t IBarkana fc Loebl 
19991: IMachacek et al.l 120011: iWise fe Abell 120071: 
O'Shea fe Normanl 120081 : iBovill fc Ricottil I2008D . We 
will therefore consider variant models (Model 3B) that 
eliminate stellar mass in pre-reionization halos below 
a critical threshold Krit.r ~ 10 km s -1 . 8 In Model 3B, 
halos with V c i TC (z TC i) < V CI it t i have stellar mass 



M* = 



(l + 0.26(K rit /Kirc(^at)) 3 ) 3 



(6) 



while halos with V^i rc (z r oi) > V cr it,r have mass given by 
equation ([5]). 

To determine very roughly the plausible range of val- 
ues for the s tellar mass fraction /* , we can refer to 
the results of IStrigari et alj (|2007| ). who derived M(< 
rtMa \)/L) = 30 — 8OOMm /iy0 for the classical dwarfs, 
and iSimon fc G cha (2003), who measured velocity dis- 
persions for SDSS dwarfs and inferred total mass-to- 
light ratios of 140 — 1800M Q /L Q . For a stellar mass- 
to-light ratio M^/Ly = 1M©/Lq, we infer plausible val- 
ues of /* ~ 10~ 4 — 10~ 2 , though these are very uncer- 
tain because all the dynamical mass-to-light ratio de- 
terminations suffer from the fact that the stars in lu- 
minous bodies of the dSphs probe only the inner parts 

8 We will refer to these as models with "pre-reionization sup- 
pression," but this simply means that halos with V c irc(zrei) below 
a critical threshold form stars with very low efficiency (too low 
to produce observable satellites), most likely because of inefficient 
cooling rather than active feedback. 



TABLE 1 
List of models used 



Model Name 



Present-Epoch Stellar Mass 



1A 
IB 



3A 



3B 



M, = /* x M sat 
Af* = /* X min((M sat /Af ) Q , 1) x M sat 

/* X M sat if V c i rc (z S at) > Krit 



M,, 



ft X M re j if V circ (Zsat) < Vcrit 



/»x(M 3at -M rci ) 



+ /. x M rci 



(1+0.26 (V crit /V circ ) J )J 

same as 3A for halos with V c irc(z r ci) > K;rit,i 
for halos with V circ (z roi ) < V crit>1 

M — /» x M Bat 

(1+0.26 (V- crit /V circ ) ii ) a 



of the dark matter potential wells. Another line of ar- 
gument comes from matching the mean space density 
of dark matter h a los to that of observed field dwarfs: 
iTinker fc Comovl (pfJQcl find /* « 10^ 3 6 at absolute 
magnitude M r w —10. In the rest of the paper, we will 
frequently refer to the stellar mass fraction normalized 
by the universal baryon fraction: 



R = 



/* 



f2(,/O r , 



6.25/* 



(7) 



Note that /* and refer to stellar fractions in 
halos where the efficiency is not suppressed, i.e., 
Vdrcf^sat) > Vcrit- We will frequently refer to the quan- 
tity (M*/M sat ) x (fi m /fi{,) as the "star formation effi- 
ciency," by which we mean the efficiency with which the 
halo converted the baryons available to it at z sat (for a 
universal baryon fraction) into stars observable at z = 0. 

3.2. Detectability and observable properties for the 
simulated satellites 

Color-magnitude diagrams for the faint dwarf 
spheroidal galaxies in the Milky Way halo show that 
the stellar populations are predominantly 'old' (older 
than several Gyrs) and metal poor ([Fc/H]< —1). To 
convert stellar masses to luminosities, we assume that 
all of our model dwarfs have a stellar mass-to-light 
ratio M*/iy w lM^/L^ appropriate to an old, metal 
poor population (jBruzual fc Charlotl l2003: M artin et al.l 
2008). The light of the lowest luminosity dwarfs can be 
dominated by a handful of bright stars and thus subject 
to stochastic variations. We ignore this complication; 
our "luminosities" are simply scaled stellar masses: 
Lv/Lq = M*/Mq. This seems appropriate, since the 
luminosities of the dwarfs galaxies are usually measured 
either by integrating over the luminosity function of old 
stellar population matche d to the observed luminosity 
function of stars in dwarfs (|Belokurov et al.l 1200 6:) or by 
averaging over possible stoch astic variations of galaxy 
luminosity (jMartin et al.l 120081) . 

The detectability of a faint stellar MW satellite galaxy 
in an SDSS-like search depends on its lu minosity and its 
distance from t he Sun, as quant ified by IKoposov et al.l 
( 2008) (see also lWalsh et al.ll2008D. O n the basis of these 
results (Figure 12 of IKoposov et alj |2008) we model the 
detectability of each simulated satellite as a binary deci- 
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Fig. 3. — Predicted stellar mass functions of all satellites within the MW's virial radius (280 kpc), for a variety of models. Left panel: 
The solid, dotted, and dashed lines represent, respectively, Model 1A with F* = 10~ 3 and Model IB with (F,,M ,a) = (10~ 3 , 10 10 Mq, 1) 
and (10 -3 , 10 10 M , 2). Middle panel: The two curves show predictions of Model 2, with F-* = 10 3 , z ro i = 11, and V C rit = 40km s 
(solid), and V cr it = 20km s — 1 (dashed). Right panel: Thin solid, dashed, and dotted lines represent Model 3A with (F* , V cr it , 2 r ci ) = 
(10~ 3 ,40km s _1 , 11), (10~ 3 , 30km s _1 , 11), and (10~ 3 ,40km s _1 ,8), respectively. The thick solid curve shows model 3B with F* = 10~ 3 , 
V ct n = 40km s — x , z Ic i = 11, and V cr i t r = 10km s — 1 . All curves reflect the average of six realizations of MW halos. These are the predicted 
complete satellite (stellar) mass functions, with no radial or sky coverage selection effects. 



sion using the criterion 

Iog lo (U /lkpc) < 1.1 - 0.228M V 



(8) 



an NFW potential with an isotropic velocity disp ersion. 
Then we can use the Jeans equation (jJean s 1919) to de- 
rive the velocity dispersion profile of stars: 



Our simulations provide the current Galactocentric dis- 
tance and orbital apocenter and pericenter for each sub- 
halo, but not the orientation of the orbit. We therefore 
assign the heliocentric distance of the satellites 



d(y(r) a 2 (r)) 



GM(r) 



= 0, 



(10) 



Dpi - a/8.5 



D 2 GC - 2 x 8.5 x Dqc cos(</>) 



(9) 



where Dqc is the Galactocentric distance (in kpc) from 
the simulations and cos(</>) is a random variable uni- 
formly distributed between —1 and 1 (0 is the angle 
between radial vectors from the GC to the Sun and 
to the sub-halo). This method assumes that the satel- 
lite orbits are isotropic ally distributed across the sky 
(see iTollerud et al~l l2008. for discussion o f the validity of 
this a pproximation). As expected from Koposov et al. 
(2008), accounting for the detectability of satellites 
causes the 'observable' population to differ strongly from 
the 'simulated' one; only the brightest satellites are ob- 
servable throughout the virial volume. 

Not surprisingly, the iKoposov et all ([2008) analysis 
also reveals a surface brightness threshold for dwarf de- 
tection, which is approximately 30 mag arcsec~ 2 with lit- 
tle dependence on distance. We assume that any model 
dwarf that passes the luminosity threshold also passes the 
surface brightness threshold. Many recent SDSS satellite 
discoveries do lie near that survey's surface brightness 
limit; this assumption can therefore only be tested with 
the next generation of sky surveys. We discuss implica- 
tions of this assumption in SJS] 

With a model that assigns stellar luminosities to each 
satellite halo, we can predict the expected stellar veloc- 
ity dispersions for comparison with th os e measured for 
MW satellit es bv I Walker et all J 200l. ISimon fc Penal 
(|2007f) . and iMartin et all (|2007D . This can be done 
straightforwardly if we assume that the stars are test 
particles — an assumption supported by the observed 
(M/L) dyn (< R cS ) > (Af/L).(< R cS ) — orbiting in 



where v is the density distribution of stars (sec 
IStrigari et all I200"7l for more detailed treatment). Here 
we assume that the density of stars follows a P rum- 
mer profile v oc [1 + (r/r p ) 2 ]~ 2 (jPlummerl fl9lTl) . which 
seems to fit observed density profiles reas onably well 
(| Wilkinson et al.ll2002t[Belokurov et al.ll2007f ). The mass 
profile M(r) used here is computed based on the virial 
radii and concentrations at the redshift z sa t of sub-halo 
accretion. Whil e the outer parts of th e sub-halos are 
tidally stripped, iPenarrubia et al] (|2008h show that the 
stars and the inner part of the dark matter sub-halo are 
stripped only at a very late stage, when the sub-halo is 
close to complete disruption. They also show that the 
velocity dispersion in sub-halos is a function of the total 
dark matter mass remaining bound inside the luminous 
body and therefore remains nearly constant until this 
late stage. 

After numerically solving the Jeans Equation, we 
compute the expected light-weighted velocity dispersion 
within the optical radius as 



J v(r)a(r) dx dy dz 
J v(r) dx dy dz ' 



(11) 



where the integration is done over a cylinder within a 
radius, R = x 2 + y 2 equal to the Plummer radius of 
the galaxy; the integral extends over ±oo in z. The stel- 
lar velocity dispersion depends on the radial extent of 
the stellar tracers, which cannot be predicted within ou r 
simple modeling context (see also iBenson et al.ll2002f) . 
We therefore use the observed properties of the faint 
Milky Way satellite s to choose stellar radii, based on 
IMartin et all (|2008l ). Specifically, we adopt Plummer 
radii r p — 150 pc for My < —5, and for fainter dwarfs 
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we adopt a linear relation between log r p and My with 
r v rising from 20 pc at My = to 150 pc at My = — 5. 

The additional important component of the detectabil- 
ity is the tidal disruption of the satellite galaxies. Al- 
though our semi-analytic model of dark matter sub-halo 
evolution properly accounts for the tidal disruption of 
sub-halos, it does not allow for the possibility that stars 
have been dispersed in a tidal stream while a small core 
of the sub-halo survives. Here we simply classify a sub- 
halo as unobservable if its current tidal radius is less than 
the expected Plummer radius of the stellar body, which 
would imply substantial tidal disruption of the stellar 
component. We also presume that a satellite is unob- 
servable if its host sub-halo has lost more than 99% of 
its original mass to tidal stripping. 

4. RESULTS 

4.1. Stellar mass function of the full satellite 
populations 

Figure [3] shows the predicted distribution of the stellar 
masses of satellites within R v iriai = 280 kpc, assuming 
47T sky coverage and complete satellite detectability. In 
the left panel, the solid curve shows Model 1A with a 
constant F* = 10~ 3 , making the stellar mass function a 
scaled version of the dark matter sub-halo mass function. 
Introducing mass-dependent suppression, Model IB with 
a = 1 (dashed) and a — 2 (dotted) lowers the low mass 
end of the stellar mass function as expected. Since this 
model also adopts F» = 10 -3 = const, above M sat = 
Mo = 1O 1O M0, the high mass end of the mass function 
is unchanged. 

The middle panel of Figure [3] shows Model 2, with 
post-reionization suppression of star formation in halos 
below a sharp circular velocity threshold, either Vcrit = 
40km s _1 (solid) or V^rit = 20km s^ 1 (dashed), where 
we have adopted = 10~ 3 and a reionization red- 
shift z IC i = 11. The resulting stellar mass functions for 
the satellite galaxies are strongly bimodal, with the low 
mass portion corresponding to dwarfs in which all stars 
formed before reionization and the high mass portion cor- 
responding to halos that exceeded the critical velocity 
threshold before becoming satellites, V^i rc (z sat ) > Vcrit- 
The low mass portion is just a scaled version of the sub- 
halo mass function at z = 2 re i- Above M* ~ 10 
the host halos are all massive enough to have star for- 
mation after z Te i, and the mass function is the same as 
that of Model 1. If the velocity threshold is lowered to 
V clit = 20 km s"\ the high mass peak in the distribution 
of satellite stellar masses extends to lower values before 
photo-ionization suppression cuts it off. 

The bimodal appearance of the middle panel of Fig- 
ure [3] is a direct consequence of the sharp V c i rc thresh- 
old for photo-ionization suppression. The right hand 
panel shows predictio ns for several v ariants of Model 
3 A and 3B, with the iGnedinl (|200l formula (Eq. Op 
used to describe photo-ionization suppression. With this 
smooth suppression, the "pre-reionization" and "post- 
reionization" portions of the mass function join to make 
a smooth overall mass function. The low mass end of the 
mass function is now a mix of satellites that formed their 
stars before reionization and satellites with V c i IC (z sat ) < 
Vcrit whose post-reionization star formation was strongly 
suppressed but not completely eliminated. Lowering the 
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Fig. 4. — Detectability of the satellite galaxies predicted by our 
fiducial model (Model 3B), as a function of their heliocentric dis- 
tance and stellar luminosity. Filled circles denote galaxies that 
can be detected with SDSS-like all-sky surveys, and empty circles 
denote those that cannot. The dashed line marks the approxi- 
mate virial radius of the MW's dark matter halo; we will com- 
pare all model predictions to the observed MW satellite population 
only within this radius. The galaxies shown were taken from one 
Monte-Carlo realization of Model 3B with (V cr it, F*, z re i, K;rit r) = 
(35 km s — 1 , 10 -3 , 11, 10km s — 1 ). The right panel shows the frac- 
tion of detectable galaxies as a function of luminosity. 



assumed reionization redshift from 



= 11 tO Z T o\ = 



1 Q- 

Conversely, if we eliminate pre-reionization SF in dwarfs 
with V" C i rc (z rc i) < V^rit.r = 10 km s 1 (thick solid line, 
Model 3B), the number of satellites with M* < 10 3 M Q 
drops by a large factor, while at higher masses the stel- 
lar mass function is unaffected. The difference between 
the thin and thick solid lines is the contribution of satel- 
lites that formed stars primarily before reionization in 
halos with V C i VC (z rc i) < 10km s _1 , for z re j = 11 and 
V cr it =40 km b -1 . 

4.2. Distribution of observed dwarf satellite 
luminosities, N{My) 

Figure |4] illustrates the impact of selection effects on 
the observable satellite population. For one realization of 
Model 3B (with parameters that yield a good match to 
observations), filled circles show satellites that would be 
detec table in an all-sky, SDSS-like survey (|Koposov et alj 
2008), and open circles show un-detectable satellites. 
The low end of the luminosity distribution, with My > 
—5, is strongly affected by the radial selection bias. 

For direct comparison with observations, we therefore 
select only those model satellites whose combination of 
luminosity and distance would make them detectable. At 
the bright end, My < —11, we assume that existing pho- 
tographic surveys are complete to Dq = 280 kpc, and we 
thus compare the total number of dwarfs across the whole 
sky to the total population of satellites within the virial 
radius in the simulation. For My > —11, we randomly 
select 1/5 of the model galaxies to mimic the 20% sky 
coverage of SDSS DR5, and we count only those satel- 
lites that would be detectable according to the criteria of 
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Fig. 5. — Model predictions for the observed satellite population, iV b s (My), including radial selection effects for the SDSS dwarfs. 
Horizontal bars show the number of currently known satellites (Table 2) in 2-magnitude bins; empty bins are plotted with an arrow. The 
SDSS and classical dwarfs are separated by the vertical line at My = — 11; note that the y-sxes for these two populations differ by a 
factor of five so that the model predictions (which incorporate a factor of 1/5 below My = —11 to account for SDSS sky coverage) are 
continuous across the boundary. Left Panel: Predictions of Model 1A, with M* <x M sa t, for three values of F*. For = 10 -4 , the green 
band shows the bin-by-bin ±lcr range of the predictions from multiple realizations; the logarithmic width of this band is similar for other 
models. Model curves have been slightly smoothed with a polynomial filter. Right Panel: Comparison of Model 1A (red curve) to Model 



IB, where the stellar mass fraction in halos with Af sa t < 10 M 







IKoposov et"aLl(|2008f) . We focus our data-model compar- 
ison on iV bs(-My"), the luminosity distribution of known 
MW satellites. We look at additional tests against stellar 
velocity dispersions, central masses, and the heliocentric 
radial distribution in § 14.31 

The luminosities, distances, and velocity dispersions 
of the observed Milky Way satellites that we use in 
all subsequent model - data compar i sons were taken 
from various auth o rs (|Mateol I1998L iMetz fe Kroupal 
I2007t iMartin et all |2Q0jj) and are compiled in Ta- 
ble [21 The sample of SDSS satellites used here con- 
si sts of those systems ab ove the 50% completeness limits 
of lKoposov et al.l (120081). We do not inc l ude two systems, 
BooI I and LeoV (jWalsh et alj 120071 ; iBelokurov et all 
l2008f ). which do not f ormally satisfy the ver y conservative 
selection limits from IKoposov et al.l (|2008f l These limits 
were chosen to avoid the issue of significant 'false posi- 
tive' detections, at the expense of leaving out 2 objects 
that deeper follow-up found to be 'real'. For the analysis 
presented here it is most important that the same selec- 
tion criteria are applied to the mock satellite observations 
and the SDSS data. As our analysis subsequently shows, 
such a small difference in sample size is smaller than 
the model halo to halo variation of number of galaxies. 
Therefore the inclusion of omission of these two objects 
does not affect our results significantly. 

Anyway, as we will sec later, the halo to halo variation 
of number of galaxies in our models is noticeable, so we 
believe that the fact that we do not include two galaxies 
should not affect our results significantly. 

The left panel of Figure [5] compares our simplest model 
(M* oc M sat , Model 1A) to the observed satellite counts, 



with a = 1 (green band) or a = 2 (blue curve). 

now including the satellite galaxy selection effects in the 
model. We randomly sample each of the six Monte Carlo 
halo simulations five times (choosing 1/5 of the faint 
satellites but always keeping the full set for My < —11), 
compute the mean model prediction as the mean of these 
30 samplings, and compute the rms dispersion among 
these 30 in each absolute magnitude bin. Despite the se- 
lection bias against low luminosity satellites, this model 
fails drastically for any choice of F*, predicting a much 
steeper luminosity function than observed. For example, 
the model with F* = 10~ 4 matches the observed counts 
near My = — 9 but predicts far too many satellites fainter 
than My = —6. Selection effects and newly discovered 
satellites have not altered th i s bas ic discr epancy, first 
empha sized by iKlvpin et al.l (|1999Pj and iMoore et al.l 
(1999). The green band shows the lrr dispersion in pre- 
dicted counts, and it is clear that statistical fluctuations 
will not resolve the discrepancy either. 

In the right panel we apply our purely empirical mod- 
ification, M*/M sat oc M a below a halo mass M sat = 
M = 1O 1O M (Model IB). With F, = 10~ 3 and a = 2, 
this model achieves reasonable agreement with the the 
observed N b s (My) over the full range > My > — 15. 
The agreement can be further improved by adjusting F« 
and Mo, so it appears that this level of mass-dependent 
suppression is approximately what is needed to explain 
the observed shape of N b s (My). Linear suppression 
(a = 1, green band) is not sufficient, predicting an excess 
of faint dwarfs when normalized to the bright dwarfs. All 
of our models fail to match the brightest bin (comprised 
of the SMC and LMC); we defer discussion of this dis- 
crepancy to the end of this Section. 
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Fig. 6. — Predictions for Model 2, in which post-reionization 
star formation is sharply suppressed below a critical velocity 
K:rit> in the same format as Figure \E\ Blue, red, green, 
and orange curves/bands show the parameter combinations 
(F*, Vcrit.Zroi) = (10- 3 ,35km s _1 ,ll), (1CT 3 , 25km s _1 , 11), 
(10 -3 , 35km s _1 , 14), and (10 -2 , 35km s _1 , 11), respectively. 
This class of models predicts a bimodal distribution of satellite 
luminosities, with the faint portion (My > —8) coming entirely 
from pre-reionization star formation. The predicted TV (My) dif- 
fers grossly from the observations. 
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Fig. 7. — Comparison of Model 2 and Model 3A, both with pa- 
rameters F, = 10 — 3 , V^rit = 35 km s , and z rc i=ll, in the same 
format as Figure [5] Switching to the continuous prescription for 
photo-ionization suppression fills in the gap between the two peaks 
of Model 2, while leaving the predictions at the highest and lowest 
luminosities unchanged. 



Figure [D shows the expected N b s (My) distributions 
for Model 2, which has a sharp V C iit threshold for the 
suppression of SF after reionization in small halos. As 
in Figure [3j the predicted N b s (My) is bimodal, with 
a bright peak corresponding to halos that exceeded V crlt 
before z sa t and a faint peak corresponding to stars formed 
before reionization. Raising the stellar fraction F* with 
other parameters fixed (orange vs. blue) shifts both peaks 
horizontally to higher My; the faint peak also increases 
in height because the brighter (though still faint) satel- 
lites can be seen over a larger fraction of the MW virial 
volume. Lowering V cr n with other parameters fixed (red 
vs. blue) has no impact on the faint peak, but the bright 
peak extends to fainter magnitudes and grows in height 
because lower mass halos can now be populated with 
stars after reionization. Raising z ro i (green vs. blue) with 
other parameters fixed has no impact on the bright peak, 
but it shifts the faint peak downwards in amplitude and 
slightly downwards in location because halos have ac- 
creted less mass by this higher redshift. While photo- 
ionization suppression reduces the discrepancy with the 
number of faint satellites seen in Model 1A, these sharp 
threshold models predict a gap between the faint and 
bright satellites that is clearly at odds with the data. 

Figure [7] compares the Mod el 2 prediction s with those 
of Model 3A, which uses the iGnedinl (|2000f ) formula to 
incorporate a smoothly increasing suppression of the stel- 
lar mass fraction in halos with V^i rc (z sat ) < V C rit- In both 
cases we use parameters = 1CP 3 , V CI it = 35km s _1 , 
z rc j = 11. Model 3A is more physically realistic than 
Model 2, with a mass-dependent suppression that is cali- 
brated on numerical simulations (and is approximately 
consistent with three independent numerical studies). 
Galaxies formed in halos with V c i IC (z sa t) < V cr it now 
fill the gap that was present in Model 2, producing 
a luminosity distribution that rises continuously from 
My = — 14 down to My = —2, before radial selection 
effects finally cut it off. With these parameter choices, 
pre-reionization dwarfs dominate the counts (and exceed 
the observations) for My < —4, but suppressed post- 
reionization dwarfs dominate the counts at all brighter 
magnitudes. 

Since Model 3 is both more physically realistic and 
more empirically successful than Models 1 and 2, we 
focus on it for the remainder of the paper, including 
Model 3B in which pre-reionization star formation is sup- 
pressed below a circular velocity threshold. Figure[5]sys- 
tematically explores the impact of parameter variations 
in Models 3 A and 3B. In the first three panels, the green 
band shows the Model 3A predictions for a fiducial set 
of parameter choices, = 10~ 3 , T4rit = 35 km s~ , and 
z Te i = 11. Changing F* (top left) shifts the predicted 
distribution horizontally to higher or lower luminosities, 
with some change in shape at the faint end because of the 
luminosity dependence of radial selection effects. Chang- 
ing V^ r it alters the predicted counts at intermediate lu- 
minosities, —4 > My > — 11, while having little effect at 
the faint end (where pre-reionization dwarfs dominate) or 
at the bright end (where most galaxies exceed the high- 
est threshold considered here). Changing z rc i alters the 
height of the pre-reionization peak at faint luminosities 
but has minimal impact for My < — 7. 

With our fiducial parameter choices, Model 3A sub- 
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Fig. 8. — Predicted iVobs(-Wv) f° r Models 3A and 3B with a variety of parameter choices, in the same format as Figure [5] In the 
first three panels, green bands show Model 3A predictions for a reference parameter set F* = 10 — 3 , V C rit = 35km S , z r ci = 11. Red 
and blue curves show the impact of changing the stellar mass fraction to F t = 10~ 2 or 10 — 4 (top left), the critical velocity threshold to 
Vcrit = 45 km s — 1 or 25 km s _1 (top right), or the reionization redshift to z rc i = 8 or 14 (lower left). The lower right panel compares 
the prediction of this reference model (now shown by the red curve) to predictions of Model 3B with a pre-reionization critical threshold 
V cr it r = 6km s _1 (green band) or 10km s — 1 (blue curve). 
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Fig. 9. — Degeneracy between F t and Vcrit for Model 3B, in 
the format of Figure [5] Blue, green, and red curves/bands show 
the parameter combinations (-F*,V cr it) = (3 X 10 — 4 , 25 km s _1 ), 
(1CT 3 , 35km s _1 ), and (3 X 10 _3 ,45km s _1 ), which all yield sim- 
ilar levels of agreement with the observations. We adopt z rc i = 11 
and V cr it r = 10 km s — 1 in all cases. 

stantially over-predicts the number of satellites with 
My w —3. Raising the reionization redshift to z rci = 14 
erases this discrepancy, but this value of z re i seems im- 
plausible given the strong and rapidly evolving opacity 
of t he intergalactic m edium at z w 6 seen in quasar spec- 
tra (jFan et al.H2006f ). and it is only marginally consistent 
with the WMAP5 results. In the lower right panel, we 
return to z rc i = 11 but suppress pre-reionization star 
formation in halos with V r c i rc (z rc i) < 6 km s _1 (green) or 
10 km s _1 (blue), motivated by the inefficient gas cooling 
expected below the threshold for atomic line excitation 
(Model 3B). The V^ r it,r = 10 km s _1 model yields accept- 
able agreement with the observed number counts over 
the full range > My > -15. The Krit.r = 6 km s _1 
model still yields an excess of faint satellites; results for 
Vcrit, r = 8 km s -1 (not shown) are nearly identical to 
those for 10 km s~\ indicating that an 8 km s _1 thresh- 
old is already sufficient to essentially eliminate the contri- 
bution of pre-reionization dwarfs. This pre-reionization 
suppression appears to be critical to explaining the num- 
ber of dwarfs observed by the SDSS. 

Within Model 3B, there is strong degeneracy between 
the values of F* and Vcrit- Figure [5] shows that the pa- 
rameter combinations (F», V C rit) = (3 X 10 -3 , 45 km s _1 ), 
(10 -3 , 35 km s^ 1 ), and (3 x 10~ 4 , 25 km s" 1 ) all yield 
similar predictions and acceptable agreement with the 
observed number counts. The lower val ues of Vcrit 
are fa vored by the numer ic al st udies of iHoeft et"aTI 
(2006) and lOkamoto et~aTl (|2008[ ). For the remain- 
der of the paper we will adopt (F„, V cr it, z r ei, Vcrit, r) = 
(10~ 3 , 35 km s~ , 11, 10 km s _1 ) as the fiducial parame- 
ter values for Model 3B. 



Fig. 10. — Fraction of pre-reionization stars in observable satel- 
lites of different luminosities, as predicted by the fiducial Model 3B. 
Filled circles show /«Af(z re j)/Af*(2 = 0), the fraction of the stel- 
lar mass that formed by z Ie i, for systems that exceeded the pre- 
reionization threshold, V c i rc (z ro i) > V cr it,r = 10km s _1 . Open 
circles show f„M(z ve i)/M t (z = 0) for systems with V c irc(zroi) < 
Vcrit, n but in the context of Model 3B these systems do not form 
any stars before reionization. The curve shows the fraction of satel- 
lites that formed more than 10% of their stars before the epoch of 
reionization, in bins of luminosity. 

For this fiducial model, Figure [TU] illustrates in more 
detail the relative importance of stars formed before and 
after reionization. For systems with Vci rc (z rc i) > V cr it,r, 
filled circles show the fraction of their stars that formed 
before reionization. For systems with V c i TC (z Te i) < V cr it,r, 
open circles show the fraction of stars that would have 
formed before reionization, but because of the V^ r jt,r 
threshold these galaxies have no pre-reionization stars 
in this model. At every satellite luminosity, the aver- 
age fraction of pre-reionization stars is small, or even 
zero, but albeit for different reasons at high and low 
luminosities. The host halos for the brighter, "classi- 
cal" dwarf satellites were typically massive enough at 
z ro i to exceed Vcrit > but that initial population of stars 
was subsequently swamped by the much larger post- 
reionization population. In contrast, the halos that now 
host the very faintest known satellites (My > —4) did 
not exceed V C rit,r at z rc i and hence — in Model 3B - 
did not form any stars before z ro i. A small fraction 
of the satellites with My — 5 have large populations 
of pre-reionization stars; these are sub-halos that just 
exceeded Vc r it, r at z ro i but have low enough values of 
Vdrc(-2sat) that their post-reionization star formation was 
strongly suppressed. If the pre-reionization threshold at 
Vcrit, r were smooth rather than sharp, then some addi- 
tional fainter systems might have significant fractions of 
pre-reionization stars. However, the general conclusion 
that pre-reionization star formation should be a small 
fractional contribution at all satellite luminosities seems 
fairly robust, provided this star formation is suppressed 
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Fig. 11. — The predicted number of MW satellites per unit mag- 
nitude within 400 kpc across the whole sky averaged from 6 MC 
realizations, using the fiducial model parameters (Model 3B with 
F t = 10~ 3 , V C rit = 35 km s _1 , V C rit,r = 10km s~\ and z rci = 11) 
and assuming no observational incompleteness. The total num- 
ber of satellites with stellar luminosities brighter than My = is 
230 ± 25. Note that this Figure gives counts in 1-magnitude bins 
rather than the 2-magnitude bins used in earlier Figures. 

in halos below the atomic cooling threshold, as seems to 
be required to match the observed luminosity distribu- 
tion. 

Figure [TT] shows the complete stellar luminosity func- 
tion of MW satellites inside 400 kpc, in the absence of 
any selection effects or incompleteness, again for the fidu- 
cia l model. (We choose 400 kpc for ease of comparison 
to llbllerud et alj [2008.) In contrast to other figures, 
it shows the luminosity function for the whole sky (47r 
sr) and in terms of dN/dMy (i.e., in bins of 1 magni- 
tude). Absent selection effects, the luminosity function 
continues to rise tow ard faint magnitudes (as noted by 
IKoposov et al. I [20081 ). contrary to the almost flat lumi- 
nosity distribution of observed dwarfs. The total number 
of satellites within 400 kpc brighter than My = ex- 
pected for the fiducial Model 3B is 230± 35. This value 
is som ewhat lower than the 400 derived by I Toller ud et all 
( 2008), but since both estimates extrapolate the number 
of known dwarfs by a factor of ~ 10, we do not place 
much weight on this difference. 

None of the models shown in Figures [5]-[5] reproduce 
the brightest observed bin — i.e., they all fail to pro- 
duce satellites as bright as the SMC and the LMC. 
Our successful models have low stellar mass fractions, 
F* ~ 10 -3 , even well above the photo-ionization thresh- 
old V^ r it . The most massive sub- halos in our Monte Carlo 
realizations have typical mass M sat ~ 10 n Mo (ranging 
from 10 10 - 5 M Q to 10 n - 4 M©), with second-ranked halos 
that are 0.2 — 0.4 dex less massive. Reproducing the 
~ 10 9 M Q stellar masses of the Magellanic Clouds then 
requires much higher stellar mass fractions F* ~ 0.05. 
To reproduce the full satellite population, the efficiency 



of gas accretion and star formation must continue to rise 
with halo mass above V cr it, or at least it must be higher 
for the SMC and LMC hosts. Since the number of bright 
SMC and LMC-like objects in our model are determined 
mainly by one parameter (because these objects are 
not suppresed by the photo-ionization) , that rise of star 
formation efficiency can not be accomodated with our 
simple model without introducing additional parameters. 

4.3. Velocity dispersions, central masses, and radial 
distributions 

As discussed in WS.21 predicting stellar velocity disper- 
sions requires assumptions beyond those needed to com- 
pute A^ b s (My). In particular, we assume that the satel- 
lites' host sub-halos have NFW profiles with concentra- 
tion given by the theoretically expected mean c(M) re- 
lation at z sa t, and that subsequent dynamical evolution 
(e.g., tidal stripping) does not alter the mass distribu- 
tion of the inner parts of the sub-halo probed by the 
stars. We also take the observed stellar radii (20 — 150 
pc, see § 13.21 for details) as input rather than predicting 
them from a physical model. With these assumptions, 
the right panel of Figure [T^] shows the predicted distri- 
bution of stellar velocity dispersions for Model 3B with 
our fiducial parameter choices. The characteristic value 
and narrow spread of velocity dispersions for the newly 
discovered SDSS dwarfs arises quite naturally from these 
models, despite the large range of stellar luminosities and 
host sub-halo masses. The predicted distribution is more 
sharply peaked than the observed one, probably because 
we did not include scatter in the halo concentration-mass 
relation and did not include observational uncertainties 
in the dispersion measurements. The mean value of a» 
differs by < 20% between data and model, but we con- 
sider this small discrepancy is not worrisome, given the 
simplicity of our dynamical modeling. 

The total masses of dwarf satellites are difficult to de- 
termine observationally because of the small extent of 
the stellar distributions relative to t he expected extent o f 
the dark matter sub-halo. However. IStrigari et alj (2008) 
show that the total mass (principally dark matter) within 
a radius of 300 pc, M300, can be inferred robustly from 
observations for nearly all of the known satellites. The 
top panel of Figure [T3] c ompares the fiducial model pre- 
dictions of A/300 to the IStrigari et alj (|2008l ) measure- 
ments. The mo del (red diamonds) n aturally reproduces 
the key result of lStrigari et alJ ([2008): over an enormous 
range of luminosities, the satellites have a narrow range 
of M300, tightly concentrated around 10 7 M Q . The the- 
oretical prediction is artificially tight because we have 
not included scatter in halo concentrations, which would 
produce roughly .15 dex (rms) of scatter in M300 (see 
iMaccio et alH 2009. fi gure 1). The model predicts a weak 
trend of M300 with luminosity, which is not evident in 
the d ata (but is similar to that predicted bv lMaccio et al.l 
[2001 . 

While the M300 range of the satellites is low, the range 
of total sub-halo masses (at z = 0) is more than three 
orders of magnitude, as shown in the middle panel of 
Figure HH The trend of total mass with luminosity is 
much stronger than the trend for M300, though there is 
a large scatter in mass at fixed luminosity because of 
tidal stripping. The near constancy of M300 is a conse- 
quence of the density profiles of CDM halos: NFW ha- 
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Fig. 12. — Predictions of Model 3B with the fiducial parameters (F t , Vcriti Zrcii Vcxit r) = (10 3 , 35 km s 1 , 11, 10 km s x ) compared to 
the observed distributions of absolute magnitude (left) and stellar velocity dispersions (right). The format of the left panel is the same as 
Figure[5] The right panel shows predicted and observed velocity dispersions only for the SDSS dwarfs — i.e., those with My > —11 — 
with data taken from Simon & Gcha (2007). 



los with the theoretically predicted c(M) relation have 
only a weak dependence of M300 on total mass over the 
range ~ 10 7 - 10 10 M W that hos ts observed Milky Way 
satellites (see IMaccio et al.l 120091 fo r further discussion) 
Thus our models and the models of IMaccio et al.l (|2009D 
are able to reproduce the narrow observed range of M300 
without much difficulty (see also lLi et al.ll2008l who ex- 
amine Meoo rather than M300). We note, however, that 
if we also allow satellites to form stars with efficiency 
-F 1 * = 10 -3 before reionization (Model 3A), then the M300 
range for the lowest luminosity dwarfs, with My > —3, 
extends downwards to A/300 ~ 1O 6 5 M0 (blue circles in 
Figure IT5]) . Thus, careful dynamical measurements for 
the faintest dwarfs could in principle distinguish whether 
they arise mainly from pre-reionization star formation or 
from highly suppressed post-reionization star formation 
in more massive halos. It is noticeab l e that our model as 
well a s the models of IMaccio et al.l (|2009D and iLi et al.l 
(|2008fl predicts that M 300 or M 600 should slightly in- 
crease with galaxy luminosity contradicting the observa- 
tions, where t here is no correlatio n at all of M300 ver- 
sus luminositv ([Strigari et al.l 120081 ). The reason of this 
disagreement is yet to be understood. It either can be 
caused by some problems with the data (selection ef- 
fects or systematics in M300 measu rements) or b y some 
astrophysical effects. For example IMaccio et al.l (2009) 
eliminates the correlation of Af3oo versus luminosity by 
assuming that the inner profile of the halos with low con- 
centration (i.e. massive halos) is modified dur ing the 
process of tidal stripping (|Kazantzidis et al1l2004f) . 

Comparing the middle and upper panels shows that 
a small number of objects have M(z — 0) lower than 
M300, which is possible because we calculate M300 based 
on the sub-halo profile at accretion. The tidal radii of 
these systems are < 300 pc, but they are all faint satel- 



lites for which the stellar Plummer radii are small. While 
their true M300 values should be M(z — 0), the values 
calculated in the upper panel are probably more directly 
compa rable to the quantities estimated by lStrigari et al.l 
( 2008) , who extrapolate to 300 pc for the faintest systems 
assuming that they are not tidally truncated within this 
radius. To minimize the tidal effects one may also com- 
pute the masses within 100 pc instead of 300 pc. For our 
simulated galaxies we also derive A4ioo, which are in the 
range 1 x 10 6 — 4 x 10 6 M Q and are also consistent with 
the M100 » 1 x 10 6 — 3 x 10 6 M Q measurements from 
IStrigari et alj (|2008| ) (supplementary information). 

The bottom panel of Fig. [T3"l shows the value of Af sat as 
a function of luminosity. The relation obviously reflects 
the underlying formula used to assign stellar masses to 
the DM halos (eqn. [5J, and the scatter caused by the 
range of accretion redshifts (which affects the M sat — V^rc 
mapping) is small. Even the faintest observable dwarfs 
have M sat ~ 10 8 5 Mq, but they have star formation ef- 
ficiencies of only ~ 10~ 5 . The difference between the 
middle and bottom panels illustrates the effect of tidal 
stripping. Nearly all the spread of M(z = 0) at fixed My 
comes from different degree of tidal stripping. 

Figure [T4l compares the distribution of heliocentric dis- 
tances of the MW satellites found in the SDSS to the 
predicted distribution for My > — 11 satellites from our 
fiducial model. We show one distribution for each of the 
six Monte Carlo halo realizations. There are significant 
halo-to-halo variations in the predicted distributions, and 
the observed distribution follows the lower envelope of 
the predictions. The distance distribution is strongly in- 
fluenced by the radial selection effects (the model predic- 
tions would be very different if we did not include them) , 
but it also depends on the radial profile of sub-halos and 
the dependence of this profile on Af sat and z S at , so match- 
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Fig. 13. — Masses of the DM sub-halos within the central 300 
pc (top), their total present-day masses (middle) and their masses 
at the time of accretion into larger halos (bottom). We only 
show halos hosting observable satellites within the MW virial ra- 
dius, as a function of satellite luminosity. Red diamonds show all 
the observable galaxies from six realizations of the fiducial Model 
3B with (F», V crit ,z rci , V cri t,r) = (10- 3 ,35km s" 1 , 11, 10km s" 1 ). 
Blue filled circles show the predictions of Model 3A, which in- 
cludes pre-reionization dwarfs (or, equivalently, has V cr it, r = 0). 
Error bars show the estim ates of M300 for observed MW satellites 
from UStrigari et al.ll2008l) . Solid lines in the bottom panel show, 
from top to bottom, M*/M sat = 10" 5 , 10" 4 , and 10~ 3 . Our mod- 
els do not incorporate scatter in the concentration-mass relation; 
adding the theoretically expected scatter would add roughly 0.15 
dex of rms scatter to the M300 predictions. 
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Fig. 14. — Comparison of the model predictions for the cumu- 
lative distance distribution of the satellite galaxies with those ob- 
served in the SDSS (black line). The predictions of the Model 3B 

with (-FV.VcritjZrei, Krit.r) = (10~ 3 , 35 km s _1 , 11, 10 km s" 1 ) are 
shown as red lines. 



TABLE 2 
Satellites used for the analysis 
and parameters adopted 



Galaxy 


Mv 




D e 


Name 


mag 


km/s 


kpc 


Bootes 


-6.3 


6.6 


60 


Canes Venatici II 


-4.9 


1.6 


150 


Carina 


-9.4 


6.8 


100 


Coma 


-4.1 


4.6 


45 


Canes Venatici I 


-8.6 


7.6 


220 


Draco 


-8.75 


10.0 


80 


Fornax 


-13.2 


10.5 


138 


Hercules 


-6.6 


5.1 


130 


Leo I 


-11.5 


8.8 


250 


Leo II 


-9.6 


6.7 


205 


Leo IV 


-5.0 


3.3 


160 


LMC 


-18.6 




49 


Sagittarius 


-12.1 


11.4 


24 


Sculptor 


-11.1 


6.6 


80 


Sextans 


-9.5 


6.6 


86 


Segue 1 


-1.5 


4.3 


23 


SMC 


-17.2 




58 


Ursa Minor 


-9.0 


9.3 


66 


Ursa Major I 


-5.5 


7.6 


100 


Ursa Major II 


-4.2 


6.7 


30 


Willman I 


-2.7 


4.3 


40 



ing the observed distribution is a significant additional 
success of the model. 

5. CONCLUSIONS 

Th e satellite discoveries in the SPSS (IWillman et all 
[200l iBebkuroy et all l200l [2007t IZucker et all 120061 ; 
llrwin et al.ll2007l ; lKoposov et all 2007; W alsh et al.ll2007f) 



have transformed our understanding of the MW's dwarf 
satellite population, extending the luminosity range by 
two orders of magnitude and the implied number of sys- 
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terns by a factor of 20. Careful quantification of the 
SPSS satellite de tection efficiency (jKoposov et al.1 [20081 : 
iWalsh et al.ll2008T) allows models that specify the relation 
between dark matter sub-halos and their stellar content 
to be tested quantitatively against the observations. We 
have shown that CDM-based models incorporating pre- 
viously advocated, physically plausible mechanisms for 
suppressing the stellar content of low mass halos can re- 
produce the observed properties of the known satellite 
population, including their numbers, luminosity distri- 
bution, stellar velocity dispersions, central masses, and 
heliocentric radius distribution. However, parameters of 
these models are tightly constrained, and alternative as- 
sumptions lead to conflict with the data. In summarizing 
our results, it is useful to review both what works and 
what doesn't. 

What works is a model in which the photo-ionizing 
background suppresses gas accretion onto halos with 



Vr.irr.( Z K!i ±) < Vcrit 



35 km s" 1 (iQuinn et all 119961: 



IThoul fc Weinberg! [19961 : IBullock et all 120001 ). with 
the smooth mass-dependent suppres sion suggested 
by numeri c al simulations (eqn. l5l IGnedinI l2000b 
Hoc ft et all [20061 IQkamoto et ail I2008D . and inef- 
ficient molecular cooling (and/or stellar feedback) 
drastically reduces the efficiency of star formation in 
pre-reionization halos below the hydrogen atomic line 
coolin g t hreshold T4 r jt. r ~ 10 km s -1 (lHaiman et aT 
iBarkana fc Loebl 119991: iMachacek et al 
IWise fc Abell |2007t IQ'Shea fc Normanl 12008: 



1997 
2001 

Bovill fc RicottT T2008) . There is some degeneracy 
between this model's two main parameters, V^ r it and 
F*, as shown in Figure but with either parameter 
fixed the other is fairly well constrained (Figure [5]) . The 
other two parameters, z rc i and V^ r it, r , just need to be in 
a range that keeps pre-reionization star formation too 
low to affect the observable luminosity function. For 
the values V cr it = 25 — 35 km s _1 favored by numerical 
simulations, F* must be < 10~ 3 , so even sub-halos 
above the V CT \t threshold have star formation efficiency 
far lower than the values » 0.1 — 0.4 found for bright 
galaxies (e.g.. iPizagno et alJ l2005t iMandelbaum et al 



2008). 



2006 1 : iDutton et all 120071 : iGnedin et al.l 120071 : iXue et al 



If we assign stellar extents based on observations, and 
make the reasonable dynamical assumptions discussed 
in §3.21 then our fiducial model naturally explains the 
characteristic value and narrow spread of stellar veloc- 
ity dis persions found for SDSS dwarfs bv I Simon fc Gcha 
(|2007f ). It also explains the charac teristic value and nar - 
row range of M 300 values found bv lStrigari et alJ {2008). 
The M300 values do not depend on the assumed stellar 
extent, and their narrow range arises from the theoret- 
ically predicted structure of CDM halos, which have a 
weak dependence of M300 on total halo mass over the 
range M ha i ~ 10 8 - 10 n M Q . Thus any CDM-based 
model that prevents formation of observable dwarfs in 
halos below ~ 10 7 Af Q should qualitatively reproduce the 
Strigari et al (200 7, 2008) results (e.g.. iLi et ail 120081 : 
iMaccio et al.ll2009l ). Tempering this success, however, is 
the fact that the total z = sub-halo masses in our model 
span three orders of magnitude; some of this range is a 
consequence of tidal stripping, but the span of M sat val- 
ues is only slightly narrower. The model, in combination 



with t he radial selection biases found bv IKopqsov et all 
( 2008), also explains the observed heliocentric radius dis- 
tribution of the SDSS dwarfs, which tests the predicted 
Galactocentric radius distribution of subhalos and its de- 
pendence on mass and accretion redshift. 

Many alternative models fail badly in reproducing 
the observed luminosity distribution. Models with con- 
stant M*/M sat predict far too many faint satellites rel- 
ative to bright satellites. The SDSS discoveries and 
luminosity-dependent selection biases do not in them- 
selves resolve the "missi ng sa tellit e" discrepancy high - 
lighted bv iKlvpin et~aTI (| 19991 ) and iMoore et all (| 19991 ): 
strong mass-dependent suppression of star formation ef- 
ficiency is still required to reconcile CDM predictions 
with observations. A simple model in which M*/M sa t = 

lO~ 3 (lVft m )(M sat /lO lo M ) 2 forM sat < 10 10 M o isrea- 
sonably successful at matching the observations. This 
successful "empirical" model has a mass dependence of 
star formation efficiency roughly like that of the success- 
ful, physically motivated photo-ionization model (eqn. O 
note that M sat cx V? TC at fixed z sat ). 

Models with sharp suppression of star formation below 
the photo-ionization threshold V cr n fail at intermediate 
luminosities, My ~ —8. Pre-reionization star forma- 
tion can provide the population of faint dwarfs in such 
a model, but there is an unacceptable gap between the 
faint and bright populations (or, for parameter choices 
that fill the gap, there is an excess of dwarfs at other 
luminosities). It is striking, therefore, that the form of 
the mass-dependent photo-ionization suppression found 
in numerical simulations is just that required to match 
the shape of the observed luminosity distribution. How- 
ever, the conversion of accreted baryons to stars must 
be very inefficient for our fiducial model to work, and it 
is not obvious why this conversion efficiency should be 
mass independent. 

The most interesting of our "negative" conclusions is 
that star formation in halos before reionization must 
be extremely inefficient to avoid producing too many 
satellites in the range > My > —6. Examination 
of Figure [8] suggests that the upper limit on the frac- 
tion of halo baryons converted to stars is a few x 10~ 4 
for z rc i = 11, or 10~ 3 i f reio nization is pushed back to 
z re i = 14. Ma dau et al.l (|2008l ) have reached exactly the 
same conclusion, with a similar numerical value for the 
efficiency limit, using the Via Lactea II simulation in- 
stead of a semi-analytic method to predict the model 
sub-halo population. Suppression of star formation in 
halos below the hydrogen atomic line cooling threshold is 
physically plausible, as the metallicity is low and molec- 
ular cooling should be inefficient. For agreement with 
Nobs(Mv), we require pre-reionization suppression in ha- 
los with Vcirc(Zrci) < V crl t,r ~ 10km S~ . 

There are several caveats to these conclusions. First, as 
discussed in H4.21 reproducing the Magellanic Clouds re- 
quires that the most massive sub-halos have M*/M sat ~ 
0.05 - 0.1, well above the F* ~ 10~ 3 of our fiducial 
model. Thus, the photo-ionization suppression described 
by equation {5| must join onto a continuing increase of 
star formation efficiency with sub- halo mass above V CI it-, 
an increase that is presumably driven by other physical 
mechanisms. Indeed, there is nothing about our results 
that necessarily picks out photo-ionization as the sup- 
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pression mechanism in low mass sub-halos, but it is a 
mechanism that comes in natu rally (one might arg ue in- 
evitably) at the desired scale (|Bullock et al ll2000{) . and 
the numerically calibrated form yields a good match to 
the observed luminosity distribution. 

In our fiducial model, even the faintest SDSS dwarfs 
form most of their stars after reionization, but they have 
K;irc(z S at) far enough below V^ r it that their star forma- 
tion is highly suppressed according to equation (O . The 
SDSS dwarfs are physically a continuum with the classi- 
cal dwarfs, and their much lower luminosities are a con- 
sequence of the highly non-linear relation between star 
formation efficiency and halo mass below V" cr it- Halos 
with Vrirc(zrei) > V CI it,r form pre-reionization stars, but 
in nearly all cases they grow large enough by z sa t that the 
post-reionization population dominates by a large factor. 
A small number of systems with My ~ —5 could have 
large fractions of pre-reionization stars, but at any lu- 
minosity such systems are rare. These conclusions are 
robust within our framework, but if we allowed for depar- 
tures from our adopted prescriptions — in particular if 
photo-ionization suppression for V^i rc <C V C mt were more 
aggressive than equation (f5]) implies and pre-reionization 
suppression weaker than we have assumed — then it 
might be possible to construct models in which many 
dwarfs with My > —6 are pre-reionization "fossils." The 
efficiency of converting halo baryons to stars in these 
systems must still be ~ 10 ~ 4 or less to avoi d prod ucing 
too many faint satell ites. [Bovill fc Ricottil (|2008l ) and 
ISalvadori et al.l (|2009D have argued that halos cooling by 
H2 before reionization naturally give rise to the physi- 
cal and chemical properties of the SDSS dwarfs. How- 
ever, eve n the low star formation efficiencies ~ 0.5% — 2% 
found by Sal vadori et al.l (1201)9; appear far too high to 
be consiste nt with the observe d number counts. On the 
other hand. lBusha et all (|2009Q propose a model in which 
post-reionization suppression of star formation is highly 
efficient (a sharp threshold) but the star formation effi- 
ciency in pre-reionization halos is strongly mass depen- 
dent, effectively spreading the low luminosity peak evi- 
dent in our Figure 6 up towards higher luminosities so 
that it fills out the entire faint end of the luminosity 
function. 

A third caveat is that we do not explain the ori- 
gin of the observed stellar extents; we just show that 
once the observed extents are adopted as inputs, then 
the observed stellar velocity dispersions emerge natu- 
rally. One possible explanation is that the baryons 
in low mass halos condense until they reach a scale 
at which the velocity dispersion is a few km s~ , and 
that this minimum dispersion provides the conditions 
necessary for star formation. We also have not at- 
tempted to explain the chemical abundance distributions 
or star formation histories of the satellites (see, e.g. , 
Orban et alj|2008t ISalvadori. Ferrara. fc Schneiderll2008l ; 
Salvadori et al.H2009f) . 

A final caveat is that we have assumed that all dwarfs 
luminous enough to be found in the SDSS also lie above 
the surface brightness thr eshold for detection, which is 
about 30 mag arcsec -2 (jKoposov et al.l I2008D . Since 
some of the known satellites approach this threshold, 
it is possible that others fall below it. A large popu- 
lation of lower surface brightness dwarfs would change 



the number counts that our model reproduces. Note 
also that a large population of pre-reionization dwarfs 
would be observationally allowed if they lie below the 
surface brightness threshold; however, even in this sce- 
nario the pre-reionization dwarfs do not account for the 
presently known satellites. Deeper large area imaging 
surveys, such as PanSTARRS, the Dark Energy Survey, 
and LSST, will show whether the MW satellite popula- 
tion includes a significant number of lower surface bright- 
ness systems. 

Our model makes several predictions that can be tested 
by these upcoming surveys or by further follow-up studies 
of known dwarfs. Deeper surveys should reveal many 
more satellites, more than 200 with My < and D Q < 
400 kpc over the full sky, with the luminosity function 
shown in Figure 111! Deep imaging of Andromeda and 
other nearby galaxies can show whether they have similar 
satellite systems, though these searches will not reach 
the extremely low luminosities that can be probed in the 
MW. Most satellites in our model have stellar extents 
that are substantially smaller than the present-day tidal 
radius of their host halo. Tidal tails and tidal disruption 
should be rare, an implication that may be challenged 
by photometric evidence on the profiles and shapes of 
the ultra-faint galaxies, which have been interpreted as 
signs of tidal distortion or disruption (e.g. Martin et 
al. 2008). Measurements of the total sub- halo masses of 
known dwarfs would provide a powerful test of the model 
predictions in Figure [T51 but the small stellar extents 
may make such measurements impossible. Our models 
predict that satellites continue to form stars down to 
z sat or below, and many observable systems should have 
z sat = 1 — 2 (see Figure [2]) . These predictions may be 
testable with detailed stellar population modeling. 

Our results gre at ly strengthe n the argument 
(Bull ock et all l2000t IBenson et all [2001 ISomervilld 
l2002t iKravtsov et al.l 120041) that photo-ionization natu- 
rally reconciles the CDM-predicted sub-halo population 
with the observed dwarf spheroidal population, thus 
sol ving the "m i ssing satel lite problem " highlighted 
by IKlypin et all (fl999l) and IMoore et all (fl999l) . The 
fiducial model presented here offers a detailed, quan- 
titative resolution of this problem in light of new, 
greatly improved observational constraints, while relying 
on previously postulated and physically reasonable 
mechanisms to suppress star formation in low mass 
halos. The MW satellites provide a fabulous laboratory 
for studying galaxy formation at the lowest mass scales, 
and much remains to be understood about gas cooling, 
star formation, feedback, and chemical enrichment in 
these systems. These issues provide challenging targets 
for numerical simulations and semi- analytic models, 
whose predictions can be tested against detailed studies 
of the dynamics and stellar populations of the known 
dwarf satellites and of the many new satellites that will 
be revealed by the next generation of sky surveys. 
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